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ABSTRACT 


The conjecture is made that the largest maximum of the normalized 


associated Legendre function P (£) lies in the interval (£,, £ ) with 

nm 1 2 


V 


^2 


|~ m(m+4) J 

n(n+l) 
^~ m(m+4) j 

n(n+l) 


where J indicates the greatest integer function. A procedure is developed 
for verifying this conjecture. An on-line algebraic manipulator, IAM, is used 
to implement the procedure and the verification is carried out for all n ^ 2m, 
for m = 1 through 6. A rigorous proof of the conjecture is not yet available. 
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I. 


INTRODUCTION 


In this report, we describe the results of an analytic study on bounds for 

the normalized associated Legendre functions P . The motivation for the 

nm 

study is to provide a rational basis for the truncation of the geopotential series 
in spherical harmonics in various orbital analyses. One would require, in addi- 
tion, bounds on the geopotential coefficients, each of which is defined as an 
integral involving the product of the Earth's density function and the correspond- 
ing spherical harmonic^. Thus, bounds on the associated Legendre functions 
may be expected to enter the truncation problem in two ways: directly as the 
functions occur explicitly in the geopotential series and indirectly as they 
occur in the geopotential coefficients. 


The goal of the study is to find a realistic upper bound for the normalized 

function P (cos 0) over the range 0 to n of the polar angle 0. The normal- 
nm 

ization factor used in this study is 


A =7(2-6 )(2n+l) , 

nm y' mo (n+m)! 

where 6 is the Kronecker 6 which vanishes for m ^ 0 and is unity for 
mo 

m = 0. The unnormalized Legendre polynomials P^feos 0) (m = 0) are known 

to be bounded between ± 1 and take on their extreme values at 0 = 0 and it . 

The normalized Legendre polynomials are thus bounded by ± «/ 2n+l. Realistic 

bounds of P for m ^ 1 appear not to be known for the full range of variation 
nm 

of 0 from 0 to ir. 


( 1 . 1 ) 


Making use of the well known integral for the unnormalized functions, 

J [p (x)] 2 dx = r^r < 2[LUB |p \V 

J - l nm J 2n+l in-m)! L 1 nm J 


-1 


a. 2) 


we obtain that for the normalized function 


1 



LUB |P (4)| n =72 

nm V mo 


(1.3) 


for m > 0. (It will be assumed that m > 0 throughout the rest of the report). 
This is a rather weak result. Another weak result may be obtained from consid- 
eration of the integral representation of the unnormalized function: 


P (cos0) = e 
nm 


im n/2 



[cos 0 + isin 0cos<p) n cosm«pd<p 


[ 2 ] 

(Hobson , p. 98, Eq. (20)). Noting that the absolute value of the integrand has 
LUB of 1, we obtain for the normalized function 


P (cos 0)1^ 
nm nm 


fo+m)J- . 

n! • n 


1 T 


-m)! (h+m)l 

> 2 

for all 0. This bound must be conservative since the fact that the integrand is 
negative over roughly half of the interval of integration is not taken into account. 

[3] 

Using an asymptotic formula from Jahnke-Emde , (p. 117), again for 
unnormalized functions, 

P (cos 0) = (~n) m y~— 2 - a" sinf(n + ^-)0+~+m^ 1 
nm V ntrsin 0 L' 2 4 2 J 

valid for 

n»m, e<0<77-e, 0<e< it/ 6 


_ J 2(2n+l)(n 


(n! 


we may conclude that in the interval (€ , n- e) the normalized functions satisfy 



This LUB, of course, becomes unbounded at the endpoints of the interval (0, TT). 

Now we know that P contains a factor sin m 0 which vanishes at 0 = 0 and 

nm 


(1.4) 

(1.5) 

( 1 . 6 ) 

a- 7) 

( 1 . 8 ) 
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0 = 77, so that the bounds for 0 near either e or 77 - e are likely to be conser- 

[ 2 ] 

vative. Hobson , (p. 302) derives and discusses a more elaborate asymptotic 
expression in which sin 0 appears in the denominator. This fact suggests that 
a realistic LUB corresponds to a value of 0 in the neighborhood of the endpoints 
of the interval (0, tt) just where the bound (1.8) is weakest. 

The associated Legendre function P (cos 0) is proportional to the 

nm 

hypergeometric function 

2 

F(m-n, n+2m+l, m+1; sin 0) 

[4] 

and it is known (see Lanczos , p. 453 and p. 369) that the hypergeometries 
functions 

1-x 


F(-r, r+2y- 1, y; 


) 


are oscillatory with the amplitude increasing as x values from 0 to 1. Hence 
if we make the identification 

r =n-m, y - m+1, x = cos 0 

we can conclude that the largest maximum of |p (cos 0) | and hence the LUB 


nm 


is 


that peak of |p closest to 0=0 (or TT, since P is either symmetric 


nm 


nm 


or antisymmetric about 0= 77/2, according as n-m is even or odd). 

In the early stages of this investigation, we examined with some care 
[3] 

Jahnke-Emde's plots of their normalized associated Legendre functions 

(pp. 112-113). Their normalization factor is just twice our A . From these 

nm 

plots it appears that for given n, | P^ | has a larger LUB than P ^^ for 
m = 2,3, ~ n, an observation which contrasts with bounds such as (1.5) and (1.8) 
which increase as m increases. 

A powerful on-line algebraic manipulator, IAM, was made available by 
Applied Data Research, Inc* for this investigation. A great many experiments 


♦Princeton, N.J. 


a. 9) 


( 1 . 10 ) 


( 1 . 11 ) 
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were tried, using IAM, and the method finally adopted was to make a systematic 

search for the smallest 0 for which the derivative of P vanishes for various 

nm 

values of n and m, with the hope of ultimately completing the analysis by a 

proof by induction. Such a value of 0 corresponds, as noted earlier, to the 

largest maximum of I P I , and the final step would be to estimate this maxi- 

nm 

mum value, again using IAM. The details of the procedure developed, the list- 
ings of the IAM programs used^nd their output appear in the following sections. 

The principal result of the study is the conjecture that the largest maxi- 
mum of P (cos 0) occurs for 
nm 


1 

" m(m+4) -1 

]* < q| 20 1 

" m(m+4) 

n(n+l) 

L 4 

J 2 n(n+l) 

.4 J 


for all m ^ 2 and all n sufficiently large - indications are that sufficiently 
large is not a stringent limitation. The conjecture has been verified for m = 1, 
2, 3, 4, 5, 6; a rigorous proof is, so far, elusive. 

The verification of this conjecture for m = 2 , 3, 4, 5, 6 involves a rather 
intricate line of reasoning. In Section II, we introduce some notations, discuss 
some pertinent properties of the associated Legendre functions, and prove the 
three theorems basic to the verification procedure. In Section III, we treat 
P as a special case, partly because it can be more completely analyzed than 
the general case, partly because some of the formalism developed for m^2 
is not well suited to the case m = 1, and partly because its relative simplicity 
can serve as an introduction to the more complex cases. Thus both Sections 
II and III may be regarded as background for Section IV. In Section IV, we des- 
cribe the IAM program used in the verification of the conjecture (1.12), and 
present the output for m = 2. The quantity of output increases rapidly with m 
so the output for m = 3, 4, 5, and 6 are contained in Appendix B. (Appendix A 
contains the IAM program and its output for the cumbersome proof of one of 
the results of Section II). In Section estimates for bounds on the largest 


* 



indicates the "greatest integer" function. 


( 1 . 12 ) 


4 



maximum of are calculated for m = 2, 3, 4, 5, 6. Once more, output from 

an IAM program is utilized. The program is straightforward and the output 

fairly lengthy; both appear in Appendix C. Finally, some comparisons are 

[5] 

made of these bounds with double precision Calculations of the largest maxi- 
mum of P , for selected values of n and m. 
nm 
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II. 


SOME USEFUL PROPERTIES OF THE ASSOCIATED LEGENDRE FUNCTIONS 


We have seen in the previous section that the bounds on the Legendre poly- 
nomials P are known. In our analysis of bounds for the associated Legendre 
no 

functions, we therefore omit the case m = 0. For all m > 0 , P (cos 0) and 

nm 

P (cos ft) vanish and the largest of the maxima of P is that one closest to 
nm nm 

9=0; this largest maximum is, of course, the LUB Ip . It is convenient to 
° 1 nm 

introduce 


t • 2 9 

? = su, - 


as the independent variable, in terms of which the normalized associated Legendre 
functions are 


P 

nm 


(£)=b (S(i-4)) m/2 

nm 


n-m 


I 


c k (-*> k 


where the coefficients C^ are given by 

_ (n-m)! (n+m+k)!m! 

k k! (m+k)! (n-m-k)! (n+m)! 


and satisfy the recursion relationship 

n 

k+1 _ x-(m+k)(m+k+l) 
C^ (k+l)(m+k+l) 

with 

C Q = 1, x = n(n+l) 


The coefficient B is 
nm 


_ = Jn+mj| __ = 1 / 2(2 

nm m! (n-m)! A nm m V ( 1} m-m)! 


(n-m)! 


We denote the polynomial factor of P by P *: 

nm nm 

n-m 


P *=Y C,(-£) 1 

nm <-> k 

k=0 


( 2 . 1 ) 


( 2 . 2 ) 


(2.3) 


(2.4) 

(2.5) 


( 2 . 6 ) 


(2.7) 
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This factor defines n-m real zeroes for P and also determines the sign of 

mn 

P for any specified value of In addition to these zeroes, the function 
dm 

P ( 4 ) possesses zeroes at 4 = 0 and 4 = 1. 
nm . 

It is rather easily shown by direct differentiation of Eq. (2.2) that 


with 


where 


for 


p 'C4) = B (4(-4)) 
nm nm 


m/2-1 - 


P * (S) 
nm 


( 2 . 8 ) 


B 

nm 


m 

2 


B 

nm 


n-m+1 


P * = 
nm 


I 


0 k «> k 


k=0 


(2.9) 



- _ (n-m)i (n+m+k-1)? (m-1)! [(m+2k)y-m(m+k)(m+k-l)3 
k k! (m+k)! (n-m-k+1) ! (n+m)! 


( 2 . 10 ) 


1 ^ k £ n-m+1 


The polynomial part P * of P * defines n-m+1 zeroes for P ’, and also 

nm nm nm 

determines its size for a given value of 4 • We note that for m=l, P 1 becomes 

nm 

infinite at 4 = 0 and 4 = 1 ; for m > 1 these points are zeroes as for P 

nm 

Now we are interested in finding the smallest zero of P * and then 

nm 

estimating the value of P at this point. Since P * and P * are of 

nm nm nm 

degree n-m and n-m+1, respectively, a complete analysis can be carried out 

analytically for n-m+1 ^ 4. In fact, a complete analysis can be carried out 

analytically for n-m+1 ^ 8, since use of the expressions for P and P f in 

nm nm 

terms of cos 0 would lead to the analysis of polynomials of the fourth degree (or 
2 

less) in cos 0 . In the present study, therefore, we could assume that we were 
seeking information in the neighborhood of the smallest zero of a polynomial 


7 



containing a fairly large number of zeroes (> 8) in the interval (0,1). We would 
thus presumably be concerned with relatively small values of the argument and 
might reasonably hope to gain some information from examination of the first 
few terms, rather than the entire polynomial. 


With such an approach in mind, we introduce some useful notation and 
set 

Pr k«> = c k-i« k ' l - c k« k 

Pr k «> = ‘W k ' 1 - c 'k« k 

L 

S L *>'1 V<-*> k 

k=0 

*L«>-Z V <*? 

k=0 

The functions Pr and Pr are pairs of terms from the polynomial parts P nm * 

and P *, respectively; S and S are partial sums for these two polynomials, 
nm ^ 

Using this notation the following properties of P * and P * maybe 

nm nm 

established: 

Property 1: The single non-vanishing zeroes and of Pr^ and 
Pr, , respectively, are increasing functions of k. Since 

K 

? k - C k-! /C k 


^k - C 'k-A 


this statement is proved if we can show that 




C k ~ C k-l C k-KL 
C k C k + l 


> 0 


( 2 . 11 ) 


( 2 . 12 ) 
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> 0 


(2.13) 


c k ' c k-i c k+i 


c k c k + i 


The proof is carried out by factoring, and discarding, the numerous common 
(positive) factors from the numerators, discarding the (positive) denominators, 
and evaluating what is left subject to the condition that 


n-m ^k+l for Pr 
n-m+1 s k+1 for Pr 

which is necessary in order for these pairs to appear at all. The calculation for 
Pr is quite easy; that for Pr is not, and was carried out using IAM. Both are 
given in Appendix A. 

Property 2: 

S T (4) >0 for 4>£ n-m s L 
L 1 

S T (4) > 0 for 4 < l , n-m+1 ^ L 
L 1 

This property follows immediately from Property 1, since all these partial 
sums are unity for 4 = 0 and cannot change sign until after the first pair be- 
comes negative. 

Property 3: If S OT ( 4 ) < 0 for some L and some 4 < C OT , then 
Z Li Z L/+Z 

P *(£) <0. Similarly if S , (4) < 0 for some L and some 4 < C OT n > 
nm 2 Li 2, Li+2 

then P nm *(l) < 0. We prove only the first statement; only minor modifications 
are required for the proof of the second. 

m- L 

= S 2L<«> - I Pr 2L + 2k<«» - C n-m« n " m 

k=l 

where £ J indicates the greatest integer function; and n-m ^ 2L. The partial 

sum S terminates with an even power of 4> and is followed by pairs of 
2 L 


(2.14) 


(2.15) 


(2.16) 
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even index which appear preceded by minus signs in P *. If n-m is even, there 
will be an integral number of such pairs; if it is odd, there will be one term left 


over which is an odd power of £ and hence has a minus sign. Now since £ < £ 


2L+2 


Pr /A = c 
2L+2k iq; 2L+2k 


2L+2k-l 

C 2L+2k 


- 0 *' 


(2L+2k-l) 


r ts >(2L+2k-li 

~ C 2L+2k ( ^2L+2k" >0 (2.17) 

Hence every term on the right hand side of Eq. (2.15) is negative and P *(£) < 0 

nm 

as claimed. 

Property 4: If S 2L+1 (£) >0 and £ < 'C 2Lf3> then p nm *(4) > 0 and 
correspondingly for P *(£). The P ro °fs are minor variations on those for 
Property 3. 


In applying these properties, we must always keep in mind that for p nni *(£)» 

n-m, the highest power of £ appearing in the polynomial, must not be less than 

the highest power of £ occurring in a partial sum, or pair, explicitly required in 

the analysis. For P *(£), n-m+1 must not be less than the highest power used. 

nm 

We are now in a position to prove the main theorems on which we base our 
calculation of bounds for the associated Legendre functions: 

Theorem I: If for a given m, there exist L, L, £ £ satisfying the 

JL ci 

following conditions: 


£ <P - 

s 2L+3 


(2.18) 


t <Ji±L 

s 2 2L+2 


^2L+2 


< C - 
s 2L+2 


(2.19) 


£ < r 

9 2 s, 2L+2 


( 2 . 20 ) 



( 2 . 22 ) 


S 2L + 1<«1> >0 


S 2L<«2> <0 


S -(£) has exactly one zero to the 
left of i 2 

S 2 -'(4)<OforO<^<4 2 

then £ and £ bracket the largest maximum of P (£) for all n such that 
1 2 nm 

n > m+2L + 2 

and this largest maximum is the only extremum of P contained in the interval 

nm 

<«!• «,>• 

Proof: The conditions (2.22) and (2.23) guarantee at least one zero in the 
interval (£ , £ ) and this zero must therefore be the one specified by condition 

X Li 

(2.24). Conditions (2.18), (2.19), (2. 22) and (2. 23), together with Properties 3 

and 4 of P *, imply that 
nm 

p nm 'l) >0 andP »(£ )<0 
nm 1 nm 2 


and hence P has at least one maximum point in the interval (£„ , £ ). 

nm 1 *2 • 

Now S -(£) differs from P *(£) by a sequence of negative Pr 

L I1IT1 

functions and, if n-m is even, the single negative term (-C ,, £ ): 

n-m+1 

5 mn*«>- 5 2L«> - I 5r 2 L + 2 k (?) 

k=l 


, „ f n-m+1 1. = >.n 

-(n-m+1 -2 — C £ 

x L 2 J n-m+1 


-m+1 


(2.23) 

(2.24) 

(2.25) 

(2.26) 


(2.27) 


(2.28) 
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Consider for a moment 


< 2 - 29 > 

This function possesses two zeroes; one at 4 = 0 an d the other at 4 = C 2 j and 
possesses a maximum at 



2J-1 *r 
2J C 2J 


(2.30) 


Since both £ and the fraction (2 J-1)/2J increase with J, n also increases 
with J. To match the parameters of the theorem, we identify J with L + 1, and 
we see that condition (2.19) guarantees that 4 0 lies t° the left of the maximum 
points of all Pr functions in Eq. (2.27). Thus, for all 4 in the interval (0, 4 0 )> 

Ct 

P *(4)-S — (4 ) is a decreasing function of 4- Further, condition (2. 25) says 
nm Z L ~ 

that S — (4 ) is decreasing in this interval, and hence, finally, P * is also de- 
2 L nm 

creasing. Since P *(0) is unity and P * (4 „) is negative.it follows that P '(4) 

nm nm 2 nm 

has exactly one zero in the interval 0<4 s 4o> so that P (4) has exactly 

2 nm 

one maximum, its largest, and no minima to the left of £ Conditions (2.20) 
and (2.21), together with the last statement, imply that P nTr t (4 ) has no zeroes 
in the interval 0 < 4 ^ 4 0 • Thus, not only do 4 and 4 0 bracket the largest 

/-A X 

maximum of P (4), but do so in a way that makes the following theorem for 
nm 

bounding this maximum useful: 


Theorem II: If the conditions of Theorem I are satisfied, then 

the largest maximum of P (4 ) satisfies the inequalities 

nm 

A < P (4) < B 
nm 


max 


P 

nm 

max 


where - 

A=v'2(2n+l)x (£ (1-4 )) m/2 P *(£ ) 

i i nm z 

B = V2(2n+l)x (4„(l-£ J) m/2 P * (4J 
2 2 nm 1 


(2.31) 


(2.32) 
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with 


x = n(n+l), n>m+2L + 2 


Proof: We know that the polynomial part, P *(4), of P (4) is an 

nm nm 

oscillatory function with n-m real zeroes in the interval (0,1), and hence n-m-1 

extrema interlaced with the zeroes in this interval. Since P *(0) =1, and 

nm 

Theorem I guarantees that no zeroes of P * lie to the left of £ , it follows 

nm 2 

that P *(4) is decreasing in the interval (0, £ ), and hence also in the interval 
nm 2 

(£l> £ 2 ), so that 


P : *(£J> 

nm 1 


p *(£j 

nm 2 


m/2 


Further, since the function (4(1-4)) vanishes at 4 = 0 and 4 = 1, and has its 
only maximum at 4 = 1/2, it is an increasing function in the interval , 4 0 ). so 

1 Ct 

long as 4 < 1/2. a condition implied by the conditions of Theorem I. Thus, 

2 

(4 l (1 "^l ))m/2 < { ^2 (1 ^2 ))m/2 


The inequalities (2.34) and (2.35) establish the conclusion of the theorem. 


We need one more theorem to make conditions (2.24) and (2.25) of 
Theorem I more explicit: 

Theorem III: Given (1) a polynomial of even degree 2K 

K 

R 2K (r)= I V k 
k=0 

with coefficients such that 

a^ > 0 k even 
a^ < 0 k odd 

and (2) r such that 

2K 

b = / a. r^ < 0 

0 £-> j 

j~0 


(2.33) 


(2.34) 


(2.35) 


(2.36) 


(2.37) 
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2K-1 


h-l 


a j+1 r J (j+l)<0 


j=0 

2K-2 


= Y a > 0 

2 L j+2 21 


b 2 = ) a 
j=0 

: 2K-i 


- 3 (.1 ) ! > 0 £ even 


£ a j +£ r £ij! <0 £ odd 


(2.38) 


b = a >0 
2K 2K 


then, in the interval (0, r), R (r) is decreasing and has exactly one zero to 

2, K 

the left of r. 


Proof: It is readily verified that the polynomial R (p) defined by 

2K 

R 2K (P) =R 2K (? + P) (2.39) 

has coefficients bj; that is, 

2K 

5 2k w)= Z y J ( 2 - 4o > 

j=o 

From Descartes' rule of signs, there is exactly one negative value of p for 

which R vanishes and hence exactly one value of r < r for which R (r) 

2K. 2K 

vanishes. This establishes the second conclusion of the theorem. 

The derivative of R (p) is given by: 

2K 

2K 

K- 2K ’< p >=I J b j P i - (2.41) 

j=0 


and from the hypotheses on the b's, the coefficients of R 2 k'^ a ^ erna ^ e s ^S n (b Q 

does not appear). Hence again using Descartes' rule of signs, R '(p) has no nega- 

tive zeroes and is therefore always negative: R *(0) = b <0. This implies that in 

2K. J. 
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the interval (0,r), R (r) decreases from R (0) = a A > 0 to R (r)=b <0. 

2K 2K 0 2K v) 

This completes the proof of Theorem III. 

In the application of Theorem III to Theorem I, we identify R (r) with 

2K 

S - ( 4 ) an d r with £ . It now remains to discuss how we evaluate the para- 

21j 2 

meters L, L, £ , £ , on which Theorem I depends. As mentioned in the 

X 2 

introduction, extensive experimentation with LAM led to the conjecture that 


„ [M-i] 

*1 n(n+l) 


V 


m(m+4) * 
L 4 


n(n+l) 

| J = greatest integer function 

will satisfy all the required conditions. To arrive at values of L and L, one 
simply starts with a small value of L, say 1, and checks to see if £ and £ , 

X 2 

defined by (2.42) satisfy all the conditions of Theorem I; if not ; increment L and 
check again. In all cases tested (m = 2, 3, 4, 5, 6), an L has been found and 
appears to be of the order of 2m. The same trial and error technique is then 
used to find L. What is lacking is a rigorous proof for the existence of L and 
L for £ and £ defined by Eq. (2.42) 

X 2 


(2.42) 
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III. CALCULATION OF BOUNDS FOR P 


iil- 


Using the notation x = n(n+l) and Eqs. (2.1) to (2.9), expressions for 
P (4) may be obtained as 

P ni ( £ } = V 2 (2n+l) [i - | + ( ^~ 2 > )( 3 X ; : - 1 4 2 ...] times JWT ) 

p »(4) V2(2n+i)xf~i -^^4 + ( x ~ 2 l . @ £ r . Q) g 2 1 — 1 
nl 2 ' L 2 s 2! • 3! J /r~ 


V4(i-4) 


We note at once from Property 4 of the last section that 


P ,(4) >0 for £< -~r = £, 
nl v>w ^ x-2 1 




so that we may identify 4^ of Theorem I with 


= 


1 3x-2 


For this case (nr=l), it is quite easy to find a value for 4 9 with L = 1, so 
that the results will hold for n ^ 2. A near-minimum value for 4 9 is obtained 
by seeking the largest a for which 


s 2 (4 2 )<0 with 4 2 


l 

x+a 


Substitution of 4„ m the first three terms of Eq. (3.2) for P ’(4) and clear- 
2 nl 

ing fractions yields 

2 

12 (x + a) - 6(3x - 2) (x + a) + (x-2)(5x - 6) 
=-x 2 +x(6a-4).+12(a 2 +a+l) < 0 


as the inequality to be satisfied by a for all n ^ 2, which implies x^6. Trial 
and error shows that a ^ 3/4 satisfies the condition, while a= 7/8 is too large. 


(3.1) 

(3.2) 

(3.3) 

(3.4) 

(3.5) 

(3.6) 
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We next verify that for a s y f > | : 

J 4 *2 1 

4 -i . -i-, >0 

2 1 x + a 3x - 2 (x + Q!) (3x -2) 

7 15 

provided x > 2 (a+1), which is so since x > 6 and 2(a+l) < 2(1 +— ) = — . Note that 
a. = 6 corresponds to the value of 4 0 given by the conjecture of the last section; 

Ct 

for a = 6 

t _I_I T m(m+4) ~| 

2 x x L 4 J m=l 


We are now ready to estimate bounds: The polynomial part of is de- 
creasing in the interval (0, 4 0 ) since £ < £ . Thus, in the interval (£ , £ ) 

Cl Li X X u 

P * has its maximum value at £ and its minimum value at £ . The other 
nl 1 2 

factor of P 7 4 (1-4 ) , increases in the interval (0, 1/2) and hence also in the 

nl 

interval (4 , £ ). Finally, therefore, the lower and upper bounds for the largest 

J- Li 

maximum of P are given by: 


A = V2(2n+l)x 7^(1-^) P nl *(4 2 ) 


B=72(2n+l)x 74 2 ( l-£ 2 ) P nl *(£ 1 ) 

respectively. We now expand the various factors in Eq. (3.9) in powers of l/x: 


(3x-2) 


(x+a) x 


x 

1+a 

x 


+ O 


(\) 


(3.7) 


(3.8) 


(3.9) 


(3.10) 


(3.11) 
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(3.12) 



*(£ 


1 


) = 1 


2E=2 _2_ + (x-2)(x-6) /_2_ \ 2_ 19 , .16 M) 

2 3x-2 12 \3x-2 / 27 81x \ 2/ 

x 


7 l+a 

p *(£)= — + — + o 
nl vs 2 ; 12 3x 


(i) 


Using these results, A and B become 

a = 7S^T) f + ^ 





In order to see that A and B not only bound the first maximum of P „ , 

nl 

but also provide reasonably realistic bounds, we tabulate A, B, and the largest 

[5] 

maximum of P calculated in double precision for selected values of n: 


n 

A 

P nl 

max 

B 

2 

1.51 

1.94 

2.23 

6 

2.43 

2.99 

3.59 

10 

'3.09 

3.78 

4.56 

14 

3.63 

4.44 

5.36 

18 

4.10 

5.01 

6.05 


(Only the first three digits of the double precision values are used for this 
table). The bounds, for this range of values for n, deviate some 15-20% from 
the exact value, which does not seem excessive for bounds intended to be valid 
for all n. Note, however, that the deviations of the bounds from the exact 
value appear to increase with n ; there is, at present, nothing in the theory on 
which to base an estimate of a rate of change in the deviations. 


(3.13) 


(3.14) 
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IV. LOCATING THE LARGEST MAXIMUM OF P 

run- 

in this section, we first describe the procedure developed for verifying 
the conjecture (2.42), then present the LAM program implementing the procedure, 
and then the output for m = 2. Because of its bulk, the output for m = 3, 4, 5, 6 is 
placed in Appendix C. 


The problem is to verify, for a sequence of values of m, that the largest 

maximum of P (£ ) lies at a point £ such that 

nm 


£ = —<£<£ = - 
S 1 x 9 x 


( 4 . 1 ) 


with 


x = n(n+l) 


m(m+ 




= z -1 


(4.2) 


The steps of the procedure are as follows: 


1. Specify m, and hence z and z^. 

2. Find L^ the minimum value of L for which 


S -<0; note n^ m + L 
2L 


3. Calculate the coefficients of u in the polynomial T - (u) defined by 

u 


LE 1 (u) = s ' 2 L 1 (f + u ) 


using Eqs. (2.9), for the coefficients C of S, and Theorem HI. Check to see 
if conditions (2. 28) of Theorem III are satisfied. If not, increase L^ by 1 and 


repeat. 


(4.3) 


(4.4) 


4. Verify that § 2 j- +1 (£]_) > 0 • 
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5. Verify that £ and £ satisfy the various other conditions imposed 

X Li 

by Theorem I for n ^ m + L . 

6. Verify that an L exists for which condition (2.23) is satisfied. 

The IAM program below, called WORK1, implements the first four steps 

of this procedure; the last two are done by "hand". IAM utilizes an ALGOL-like 

language. It executes PARTS of a program and/or STEPS within the parts. The 

individual command strings are identified in a "decimal" notation: PART. STEP, 

eg. 4.3, 5. 11, etc. Note that STEP 5.11 would be executed between STEPS 

labeled 5 . 1 and 5 . 2 . The commands are almost self-explanatory; T denotes 

exponentiation, is a replacement operator, FOR generates a loop, WITH 

indicates substitution. All arithmetic is carried out in integer form to full 

precision, using as many "words" as necessary to store the result. There are, 

of course, limits to the storage capability, and when the address space of the 

system is exhausted, a message to this effect is returned and calculation stops. 

WORK1 successfully carried out the first four steps of the procedure outlined 

above for m = 2,3,4, 5, and 6. Time did not permit further trials, and most 

likely storage would be exhausted for m = 7 or 8; considerable care in conserving 

storage had to be exercised to carry through m = 6, for which 14 terms of the 

P ' series were required, 
am 

With these preliminaries, here is WORK1: 


TYPE WORK 1 

1*1 : 12- I 1 - 1 5Z1-Z- 1 5 TYPE M. I I > 1 2 . Z. Z 1 
1.2:MSQ-Mt2 
1 .3 :TW0M-2*M 
1. 4 JMMC 1 ) -MSQ + M 

1.5:F0R 1-2 TO 1 2t MMC I >-MM C I — 1 > +TWOM + 2* (1-0 
1 • 6 : MSQ-2 *MSQ 

1 • 7 : FOR 1-1 TO 1 2., MMM ( I >- ( 1 + I ) *M + C I + 1 > * ( I + n 
4. 1 :O(0)-MJ FOR 1-1 TO I 1 , 0 < I > - 0 C I - 1 > * J * ( M+ I > ; & 

CPC05-QCI 1 ) ; DELETE Q 
4.2 :CP< 1 >-CPC0)/MM< 1 ) 

4.3: FOR 1-2 TO I 1 > (CP ( I >- < CP< I - 1) /MMM C I - n > * ( X-MM C I - 1 > ) « 
.•DELETE MMM ( I - 1 ) ) 
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5. 1 :MMM< 1 )-M*MMC 1 ) J MMM ( 2 > -M S0+ TWOM 
5. 1 5 : DELETE MM 

5.2 : FOR 1-3 TO I 1 , MMM < I) -MMM < I - 1 ) +TWOM 

5.3 :Y( 1 >-<M+2)*X-MMM( 1); DELETE MMM ( 1 > 

5. 4 : FOR 1-2 TO I 1 , ( Y C I) - Y C I - 1) + 2* X-MMM ( I > & 

^DELETE MMM C I ) ) 

5. 6 : FOR 1-1 TO I 1 , ( CP C I > -CP C I > *Y C I > * ( - 1) » I i & 

DELETE Y(H) 

6. l :D<0>- 1 ;D( l)-x 

6. 2 : FOR 1-2 TO I 1 , D C I >-D C 1 > *DC I - 1 ) 

6. 3 : FOR 1-0 TO I 1 , CP C I ) -D C I 1 - I > *C PC I) 

6. 4 : DELETE D 
8* 1 : $D I ST- 1 5 

8*3 : X2-U+ C I2 + M) * ( I2+M- 1 ) 

8*4:Xl-U+( I 1 + M ) * C I2 + M) 

8»5:AI 1-SUM ( 1-0 TO I 1 :CPC I ) *Z 1 ♦ I ) 

8. 6 : ODD- A I 1 WITH [X=X13jTYPE ODD 
8.65: DELETE All, ODD 

8. 7 :AI2-SUM( 1-0 TO I2:CPCI)*Ztn 
8.8:TESTC0)-AI2 WITH CX=X23iTYPE TEST ( 0 ) 

8. 9 : DELETE TEST(0),AI2 

9. 1 : FOR J-l TO I2,IPCJ)-J 

9.2 :T-SUM( I- 1 TO 1 2 : C I P C I > * Z » C I - 1 > > *C PC I > ) 

9.3:TESTCn-T WITH CX=X23 

9 . A : TYPE TEST C 1 ); DELETE TEST(l) 

9 . 5 : FOR K-2 TO CI2-l>>CFOR J-K TO I 2 » I P C J > - I P C J > * & 

CJ-K+1 )/K; T-SUM C I -K TO 1 2 : C I P C I ) * Z t C I - K 3 3 *C PC I) > ; & 

TEST ( K )- T WITH [X = X23,*TYPE TEST ( K >; DELETE TEST C K > , T ) 

9.6 : DELETE ALL VALUES 
10. 1 :C (03-CPC0)/XT I 1 

10.11 : FOR 1-1 TO 12, CC I 3-COEFFCCPC I ) , X, I 1 ) 

10. 12 :TST1 (03-SUMCI-0 TO I 2 : C C I ) * Z ♦ I > i TYPE TST1C0) 

10. 1 3 : I F TST1 (0) . GT. 0, ( I 1 -I 1 +2; DELETE TSTliTO STEP 21.1) 
10.15 :FOR J-l TO I2,IP(J>-J 

10.2: T-SUM C I- 1 TO I 2 : C I PC I > *Z r C I - 1 ) ) *C C I ) ) 

10.3 :TST1 ( 1 )-T 
1 0 . 4 : TYPE TST 1 C 1 > 

1 0 . 5 : FOR K-2 TO CI2-1>,(F0R J-K TO I 2, I PC J) - 1 PC J ) *& 

(J-K+l )/Ki T-SUM C I -K TO 1 2 : C I P C I ) *Z t C I -K ) ) *C C I) ) J & 

TST 1 ( K ) - T J TYPE TST1CK);IF C TST 1 ( K ) * ( - 1 > t K ) . LT . 0, C I 1 - 1 1 +2 i & 
DELETE TS T 1 ; TO STEP 2 l . 1 > > 

10.8:DELETE C,T,TST1 
10. 9 :DO PART 8 ) DO PART 9 
21 . 1 : 12-11 - 1 

21 . 3 : DO PART 1 J DO PART 4 

21 . 4 : DO PART 5iDO PART 6JTO STEP 10.1 
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Execution of the program is accomplished by entering the command 
DO PART 21 from the terminal; this must be preceded by entry of m, z,and 
an initial value for an odd integer II. The even integer 12 = 11-1 plays the role 
of 2L in step 2 of the procedure outlined above. Note that PART 21 executes 
PARTS 1, 4, 5, and 6 and then skips to STEP 10.1. PARTS 1, 4, 5, and 6 evaluate 
an array CP of (II +1) coefficients which are proportional to the first (II +1) 
terms of the polynomial part of P m '(4). The proportionality factor is the 
least common denominator of these coefficients, including the power of x 
which will occur on substitution of £ or £ for the argument. IAM works 
more efficiently if it is not asked to carry fractions, since its first step is to 
get the common denominator, which requires storage and is carried throughout 
the calculation. 

Referring to Eqs. (2.8), (2.9) and (4.1), we note first that 




(4.5) 



.[x-m(m+l)][x-(m+l)(m+2)j»»*[x-(m+U-2)(m+Il-l)] . [(m+2Il)x-m(m+Il-l)(m+Il)] ? 


II! m(m+l)». • (m+Il) 


and the corresponding terms in P * are these coefficients multiplied by the 

nm 

corresponding powers of z/x or z /x for evaluation at £ and £ respectively. 

1 J- " _ 

In either case, the lowest common denominator of the partial sum S , which 

... - - II 

serves also as a common denominator for S^, will be II! m(m+l)» • • (m+Il)x 

and CP (I) at the end of PART 6 of WORK1 is this factor times the corresponding 
Cj of Eq» 4.5: 
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(4.6) 


CP(0)=I1! m(m+l)- • • (m+Il) x 1 

C P(1) = -I1! (m+2)(m+3). • • (m+Il)x^[(m+2)x-m 2 (mHT)] 

CP(I1) = (-l) I:L (x-m(m+l))(x-(m+l)(m+2))--* 

(x-(m+Il-2) (m+Il -1)) ((m+2Il)x-m (m+LL-l) (m+Il)) 

PART 1 and STEPS 4.1 and 4.2 set up the factors independent of x in these 
expressions. STEP 4.3 calculates all but the power of x and the factor in- 
volving (m+2I)4. The array Y is set up in STEPS 5.3 and 5.4 to generate 
this last factor and STEP 5.6 inserts it along with the proper sign into the 
CP's. Finally, in PART 6 the array D generates the powers of x which are 
inserted into the CP's by STEP 6.3. 

Having generated expressions for the terms of the partial sums, we 
next want to find of Eq. (4.3). The search is done by successively calcu- 
lating 

12 

^12 I CP(k) ^ (4 * 7) 

k=0 

for 12 = 2, 4,6, . . . until a negative value results. In practice, this is a very 
lengthy calculation which strains the storage capacity of IAM, so a similar pre- 
liminary search is carried out using only the leading terms in x. Note from 
Eq. (4.6) that each CP has leading term x 11 and that x = n(n+l) can be arbitrarily 
large. Thus, in effect, we first seek an for very large x. This is accomplished 
in STEP 10.11 by constructing an auxiliary array C(I) from the coefficients of 
in CP(I). The rest of PART 10 calculates and prints out 

12 

TST1(0)=£ C(k)z k (4.8) 

k=0 
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incrementing 12 by 2 and repeating the calculation until a negative value 
results. At this point the test is undertaken for zeroes of this sum to the left 
of z, using Theorem III. If more than one zero is detected, 11, and hence also 
12, is incremented by 2. The sequence of tests for Theorem III requires 


TST1(0) <0 
12 

TST1(1)=^ kC(k)z k <0 
k=l 

TST 1 (2) = V C(k)z k >0 

i—k 1.2 

k=2 


TST 1(12-1) = - C(I2-l)z 12 1 +I2C(I2)z 12 <0 
TST 1(12) = C^z 12 >0 


The last of these tests does not appear in the output. Once these tests are passed 
for the leading terms in x, the value of 12 so obtained is used as the initial 
value for the corresponding sequence of tests using the complete coefficients of 
Eq. (4.6), with the following modification: 

First, note that use of 12 for P * requires 

nm 

n-m+1 £ 12 


and hence 


x ^ (m+I2-l)(m+I2) 


WORK 1 defines 

x2 = (m +12) (m +12-1) +u 


xl = (m+Il) (m+I2) + u 

with u ^ 0 to incorporate the fact that u = 0 corresponds to the minimum value 
of n(n+l) consistent with the existence of Sjg and S^, with 12 corresponding 
to 2L and II to 2L+1 of Theorem I. In STEPS 8.5 and 8.6, IAM first 


(4.9) 


(4.10) 

(4.U) 

(4.12) 
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evaluates S^(z^), labeled ODD in the output, with x = x^. If the coefficients 
of all powers of u are positive, then will be positive for all n ^ m+Il-1, 

as required by Theorem I. In STEP 8.8, S (z), labeled TEST(O) in the output, 

I & 

is evaluated with x = x . If the coefficients of all powers of u are negative, S 
will be negative for all n ^ m+I2-l, again as required for the application of 
Theorem I. In PART 9, the calculations corresponding to Eq. (4.9) are carried 
out using the full coefficients CP, x=x~, and arranging the output in powers of 
u. The results of the calculation are returned in the array TEST(J), J = 0, 1, . ..12. 
The inequalities corresponding to those of (4. 9) will be valid for all n ^ m+I2-l if, 
and only if, the coefficient of each power of u satisfies the inequality individually. 
This test is made by "human" inspection of the output. In all cases run, the 
value of 12 found by the preliminary tests led to an array TEST(J) satisfying 
the above conditions. 

At this point we display the output of WORK1 for m=2. For m=2 

z = [lip>3 

and we initialized II to 3: 

• I AM 

WELCOME TO I AM C 7232 1 ) 

♦LOAD FROM "WORK 1 " 

*M-2J I 1-3 J Z-3iD0 PART 21 

M: 2 


I 1 t 3 
12 : 2 


Z: 3 


Z1 : 2 

TST1C0): 90 
M: 2 
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II: 5 



12: 4 


Z: 3 
Z1 : 2 

TST 1 ( 0 ) : - 77490 
TSTl(l): - 93240 
TST 1(2): 34020 
TST 1 < 3 ) : - 2520 
ODD: 

52256*11* 5 + 1 1 387840*1)* 4 

+ 9928 1 9328*11 * 3 + 4328829 1 584*11* 2 + 94399 4926080*U + 8237259878400 
TEST ( 0 > : 

- 77490*11* 5 - 1 1872980*11*4 

- 721700280*11*3 - 21810222000*11*2 - 328224960000*11 - 1969349760000 
TEST ( 1 ) : 

- 93240*11*5 - 14802480*U*4 

- 930696480*11*3 - 2909 1 1 82 400* U* 2 - 453309 69 6000*11 - 282 17 54880000 
TEST < 2 > : 

34020*1) * 5 + 4906440*11*4 ~ 

+ 287083440*11* 3 + 847370 1 600*U* 2 + 125574624000*U + 744629760000 
TEST (3 > : 

- 2 520*1) * 5 - 408240*U* 4 

- 24625440*11*3 - 706708800*U*2 - 9761472000*U - 52254720000 
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We see that Eq. (4.13) led to a value of 90 for TSTl(O), so WORK1 incremented 
11 (and hence also 12) by 2 and started over. This time the array TST1 pro- 
duced the acceptable pattern as noted above. TST1(I2) must be positive, 

unless the last coefficient CP(I2) is incorrectly computed. Next ODD is calculated 
and we see that it is positive for all non-negative u. Finally, TEST(O), TEST(1), 
TEST (2) and TEST (3) are calculated, the output is inspected, and it is verified that 
for u £ 0, the sign pattern -, -, +, - is correct. 

At this point we may conclude the following: For m=2 

S 4 (f) <0 n*5 

S f-)>0 n ^ 6 

5 NX / 


Referring to Theorem I, we identify 





We must next verify that conditions (2.18) and (2.19) are satisfied. Before doing 
this we discuss three special cases, corresponding to n = m+2L^-l, n = m+2L^, 
n = m+2L^+l, for which these conditions may not be relevant. The existence of 
the partial sum S — implies that the first of the special cases represents the 

2 Li j ^ 

smallest value of n consistent with given m and L^. 

1. n = m+2L -1: S - =P * and therefore neither S - ■ nor the 

1 2 Li^ nin 2 JLi_^ * JL 

Pr’s associated with conditions (2.18) and (2.19) appear. These conditions must 
be replaced, in this case, by the condition 



(z x ) > 0 


2. n = m+2L 4 : S— = P *. In this case conditions (4.14) replace con 

1 2L^1 nm ’ ' 

ditions (2.18) and (2.19): p nm *( z ) = S jr ~ ^i-m+1 Z n_m+1 < 0 if S - (z)<0. 

nm n-m+ 2L 1 

3. n = m+2L +1: In this case 


(4.14) 


(4.15) 


( 4 . 16 ) 
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(4.17) 


(2Lj_+2) 


f nm* (Z) = S 2L, “ Pr 2L 


= s - + C - z 

2 L l +2 2^+1 2 L x +2 


Hence condition (2.19) is required, but (2.18) is not, since if Sg- >0, then 
also P * > 0. * 


nm 

Verification that the modified conditions hold for these special cases is 
easy. We proceed to the verification of conditions (2.18) and (2.19) for the 
general case for m = 2 ; n ^ m+2L^+2: 


p _ p S = 8-10 (16X-2-8-9) = 640 1-9/x 2 

SjL+3 £ (n-8)(n+9)(18x-2* 9*10) 9x (l-72/x)(l-10/x) x 

8 

C 

p 6 _ 7° 9 14x-2* 7» 8 = 441 1-8/x > 3 

aLj+2 C ? (n+8)(n-7) 16x-2-8-9 8x (1-50/x) (1-9/x) x 


(4.18) 


for all n ^ 8, which implies x ^ 72. 


We have now satisfied all conditions of Theorem I, except that on S (£ ). 

2 Li 2 

This last test is carried out by another IAM program EVAL, which is used, not 
only for this test, but also in the next section, for the calculation of bounds using 
Theorem II. Here we shall merely state that EVAL does indeed verify that L 
exists for which conditions (2.20) and (2.21) of Theorem I are satisfied. 

The output of IAM for m = 3,4,5, and 6 appear in Appendix B. Because 
of prior experimentation , was already known for these values of m, and 
to reduce the bulk of output, 11 was initialized to 2L^+1. The output exhibits all 
the desired "sign" patterns for the array TST1, ODD, and the array TEST for the 
values input for m, 11, z, z^. Verification of conditions (2.18) and (2.19) for 
the general case, or their substitutes for the special cases, is routine, though 
tedious; the details are omitted from the report. 


The table below summarizes the information given by the output of 
WORK1 (plus the additional verifications), showing 12 = 2L , £ , £ and the 

i. X z 

lower bounds on n and x = n(n+l) for which £ and £ bound the largest 

X z 

maximum of P : 
nm 
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n 

12=2 


^2 

lower 

n 

Dound for 

X 

2 

4 

2/x 

3/x 

6 

42 

3 

4 

4/x 

5/x 

n 

56 

4 

6 

7/x 

8/x 


110 

5 

8 

10/x 

ll/x 


182 

6 

12 

14/x 

15/x 

m 

342 
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V. 


BOUNDS FOR THE LARGEST MAXIMUM OF P 


nm— 

The case m = l is discussed in Section III. For m = 2, 3, 4,5,6, 4^ and 
4 given in the table at the end of the last section can be used, together with 
Theorem II (Section II) to obtain the desired bounds. For this purpose, we re- 
quire approximate values of P *( 4 ,) and P *( 4 „) for m = 2, 3, 4, 5, 6. 

nm 1 nm 2 

These values of the polynomial parts of P nm were calculated using the LAM 
program EVAL, reproduced in Appendix C. The code is fairly obvious. The 
output arrays SUMX(I) and SUMIN(I) refer to (2m)! y m times the first I 
terms of P * evaluated at 4, and 4„> respectively. The parameter y re- 
fers to n(n+l), denoted by x in the rest of the report. It will be noticed that 
for m = 2, the numerator of 4 0 » z > was inadvertently set equal to 2, instead 
of 3 in the first command. We shall therefore outline here the calculation for 
m = 2, and this will serve as an illustration of how EVAL calculates SUMX and 
SUMIN for m = 3, 4, 5, 6. 

From Eq. (2.3) 

P ♦, £v=1 _ .(n- 2 ), fa + 3) ... 

n2 3 9 2! -3-4 

from which we readily obtain 


p * \ i - x - 6 t , (x-6)(x-12) t 2 (x-6)(x-12)(x-20) .3 

Jr _o VS / A o’ o A S ~ trrjo S ••• 


n2 
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5-72 


Substituting 4 -, = 2/x an d = 3/x and arranging in inverse powers of x, we 

X Z 

obtain 

43 83 12 160 

SUMX(3, = - + — + — + — 

5x 3x 

3 21 27 108 

SUMIN (3) = T7T + ~ TT - 

w 10 lOx c 2 3 

5x x 

O 

We now observe that, while powers up to 4 were retained in calculating 
SUMX and SUMIN, Theorem I requires n ^ 6 and hence x s 42, so that only 


(5.1) 


(5.2) 


(5.3) 
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the first two terms need be retained for approximate values of P *: 

n2 

P *(£ ) ~ 
n2 ' s 2 ; 10 lOx 

n2 v *r 90 45x 

To implement Theorem II, we next calculate B in powers of l/x. From Eq. 
(2.6), it is readily verified that 

B = — \/2(2n+l)x(x-2)(x-6). . . . (x-m(m-l)) 
nm m! 

= 2x m//2 (n+l/2)(l-2/x)(l-6/x). . . (l-m(m-l)/x) 
m! 


and for m = 2 


n2 


= x J (n+l/2)(l-2/x) = Vn+l/2 x [l - - - . . „ J 

X 2x 


~Vn+l/2 x (1-1/x) 

2/2 

We also need to evaluate (£ (l-£)) at £ and £ • 

X ^ 

«2 (1 -«2 )= f( 1 -f) 

We now insert these various expressions into Eqs. (2.32) for A and B, the 

bounds on the largest maximum of P we obtain 

n2 

(i-f) 


\ x/ x ' x/ \ 90 45x 


— V n+l/2 (43-6/x) 


(5.4) 


(5.5) 


(5.6) 


(5.7) 


(5.8) 


31 



as approximate values. 

the largest maximum of 
selected values of n: 


The table below tabulates the values as well as 
P by a double precision calculation 


n2 


[ 5 ] 


max 

for 


n 

A 

P n2 

max 

B 

4 

.64 

2.16 

3.02 

5 

.94 

2.35 

3.35 

6 

1.16 

2.53 

3.64 

10 

1.77 

3.18 

4.64 

15 

2.26 

3.84 

5.64 

20 

2.65 

4.41 

6.49 


Note that we have included the values n = 4, 5, even though they are below the 

bound imposed by Theorem I. We do the same for m = 3, 4, 5, and 6 below, since 

it appears that Theorem I is conservative in this respect; in fact, some preliminary 

analysis to generalize Theorem I indicates that the bracketing of P nrn by £ ^ 

max 

and £ holds for all n ^ 2m. 

2 

The output of the arrays SUMX and SUMIN from EVAL is included in 
Appendix C, along with the listing for EVAL, for m = 3,4,5, and 6. Tables 
similar to that for m = 2, above, were obtained by the following sequence of 
steps : 

1. Replace y in SUMX and SUMIN by x. 

2. Approximate P *(£„) and P *(£„) by truncating 

run 1 nm ^2 

p v SUMDjiml . 

nm ls l/ ~ • ' m 


P 

nm 


*«a> 


SUMX(m) 
(2m)! x m 


(5.9) 
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to a reasonable number of terms in l/x, recalling the lower bound for x given 
at the end of Section IV. 

3. Expand (4 1 (1-4 )) m ^ 2 and (4 0 (l - 4 9 )) m ^ 2 > including all powers of 

XX Ct Ci 

l/x retained in 2. 


4. Expand B^^, as given by Eq. (5.5) similarly. 

5. Calculate A and B from Eq. (2.32), again in truncated form. 

6. Evaluate A and B for the desired values of n and x=n(n+l). 


The tables so constructed are: 


m=3: 


n 

A 

P 

n3 

max 

B 

6 

2.10 

2.30 

3.34 

9 

2.05 

2.72 

3.79 

12 

2.15 

3.10 

4.22 

15 

2.29 

3.44 

4.65 

18 

2.51 

3.75 

5.18 

20 

2.53 

3.95 

5.30 


m=4: 


n 

A 

P n4 

max 

B 

8 

2.25 

2.42 

3.30 

10 

2.18 

2.65 

3.48 

12 

2.28 

2.87 

3.76 

15 

2.51 

3.18 

4.18 

17 

2,69 

3.37 

4.46 

20 

2.90 

3.64 

4.84 
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1X1=5: 


m=6: 


n 

A 

P 

n5 

max 

B 

' 

10 

2.13 

2.51 

2.78 

11 

2.28 

2.61 

3.03 

13 

2.37 

2.81 

3.41 

15 

2.60 

2.99 

3.69 

18 

2.74 

3.25 

4.05 

20 

2.83 

3.42 

4.26 


n 

A 

P „6 

max 

B 

12 

1.98 

2.59 

3.26 

14 

2.28 

2.76 

4.43 

16 

2.44 

2.93 

4.23 

17 

2.51 

3.01 

4.21 

18 

2.57 

3.09 

4.22 

19 

2.63 

3.17 

4.24 

20 

2.69 

3.24 

4.28 


The principal conclusion to be drawn from these tables is that the values provided 
for £ and £ , from the conjecture of Section II, appear to provide realistic 

X ct 

bounds on the largest maximum of P 

nm 
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APPENDIX A 


PROOF OF THE INEQUALITIES (2. 13) 


Proof of the Inequalities (2.13) 

The first of these inequalities is 




> 0 


where the C's are coefficients in the polynomial part of P . For all three 

nm 

of these coefficients to appear in P , k must exceed zero. Since the C's 

nm 

are all positive 




> 0 


is an inequality equivalent to (A.l), but somewhat simpler to manipulate. From 
Eq. (2.4) we have 


and hence 


with 


Q 

k+1 _ x- (m+k) (m+k+1 ) 
C k ~ (k+l)(m+k+l) 


^k _ x- (m+k-1) (m+k) 
C k _ 1 " k (m+k-1) 


x = n(n+l) s (m+k+l)(m+k+2) 


from the. first of inequalities (2.14). Substitution of Eqs. (A. 3) and (A. 4) into 
Eq. (A. 2), and clearing fractions, yields the following inequality to be verified; 

2 r -i 

(m+l)x-2m(m+l)(m+2) + (m+1) (m+2) = (m+1) [x-(m+2)(m-l)] > 0 


From inequality (A. 5) 

x-(m+2) (m-1) ^ m^ + (2k+3)m+(k+l)(k+2)-(m^+m-2) = 2(k+l)m+(k+l)(k+2)+2 > 0 


(A.l) 

(A. 2) 

(A. 3) 

(A. 4) 
(A. 5) 

(A. 6) 
(A. 7) 


A-l 



which establishes (A. 6) and hence (A.l). 


The second inequality to be verified is 



C 



k 




>0 


where the C’s are coefficients in the polynomial part of and, as above, 

k 2 : l is required if all three are present. Also, as above, all the C's are 
positive and hence the denominator in inequality (A. 8) may be discarded. The 
case k=l is a special case for which (A. 8) reduces to 


since, according to Eq. (2.10), C Q = 1. and are given by 

- _ (m+2)x-m 2 (m+l) 

C 1 m(m+l) 


- _ (x-m(m+l))((m+4)x-m(m+l)(m+2)) 
2 2! m(m+l)(m+2) 


The verification of inequality (A. 9) is laborious, but has been carried out by 
hand since it was inadvertently omitted from the LAM program for the general 
case. The details of the proof are omitted; the key step is to make use of the 
second inequality of (2.14) 

n ^ m+k =m+l, which implies n = m+L+1 with L ^ 0 


so that we may set 

x = n(n+l) = (m+L+1) (m+L+2) 

where L is any integer ^ 0. Then the sequence of steps 

1. Substitute (A. 10) into (A. 9) 

2. Clear fractions 

3. Collect terms in ascending powers of x 


(A. 8) 

(A. 9) 
(A. 10) 

(A. 11) 
(A. 12) 



4. Substitute for x from (A. 12) 

5. Arrange in ascending powers of L 

6. Observe that the resulting expression is non-negative for non- 
negative values of L. 

For the general case, k > 1, the three coefficients involved in (A. 6) are, 
from Eq. (2.10) 

- _ (n-m)! (n+m+k-2)! (m-l)i [(m+2k-2)x-m(m+k-l)(m+k-2)] 

k-1 (k— 1) ! (m+k-1)! (n-m-k+2)! (n+m)! 

- _ (n-m)! (n+m+k-1)! (m-1)! [(m+2k)x-m (m+k) (m+k-1)] 

k k! (m+k) ! (n-m-k+1)! (n+m)i 

- _ (n-m)! (n+m+k)j (m-1)! [(m-t-2k+2)x-m(m+k)(m+k+l)] 

k+1 (k+1)! (m+k+1)! (n-m-k)! (n+m)i 

We discard the denominator and the positive factor 

(n-m)! (m-1)! 

(n+m) ! 

common to all three C's. We next clear fractions, discard the positive common 
denominator, and a number of other common positive factors in the numerator 
to obtain, as the inequality to be verified 

(k+1) (m+k+1) (n-m-k+2) (n+m+k-1) [(m+2k)x-m(m+k) (m+k-1) 

-k(m+k) (n-m-k+1 ) (n+m+k) [ (m+2k-2)x-m (m+k-1) (m+k-2) ] 

times [(m+2k+2)x-m (m+k+1) (m+k)] >0 

The proof of this inequality is very laborious and LAM was used to carry out 
most of the remaining calculations. The program written for this purpose 
follows : 


(A. 13) 


(A. 14) 


(A. 15) 


A -3 



TYPE PROOF 

30.1 :N-M+K+Jj TDI ST- 1 5 
30.2:F1-K*(M+K)*(N-M-K+1)*(N+M+K> 
30.3:F2-(M+2*K-2)*(N+1 )*N-M*(M+K-1 ) * CM + K-2) 
30. 4 : F3- (M+2*K+2)*N*(N+ 1 )-M*(M+K+l ) * ( M + K > 
30.5:F4-(K+1)*(M+K+1)*(N-M-K+2>*CN+M+K-1) 
30.6:F5-((M+2*K)*(N+l)*N-M*(M+K)*(M+K-l))t2 
30.7 :TST-F4*F5-F1 *F2*F3 

30.8:TSTJ(0)-TST WITH [J = 0UTYPE TSTJC0) 

30 . 9 : FOR 1-1 TO 6 > C TST- ( TST- TST J ( I - 1 ) ) / Jj & 
TST J ( I ) -TST WITH [J=0];TYPE TSTJ(I)i& 

DELETE TSTJ ( I - 1) ) 

30. 95: DELETE TST , TST J, N , F 1 , F2,F3> F4>F5 


The first instruction makes use of condition (2.14): 

n ^ m+k 

so that J in the program must satisfy 

J ^ 0 

The command $DIST <- 15 requires the system to expand all powers of multi- 
nomials up to and including the 15th. In terms of the parameters FI, F2, F3, 
F4, F5 generated in STEPS 30.2 to 30.6, the left hand side of (A. 15) is given by 

F4*F5-F1*F2*F3 

as indicated in STEP 30.7. This expression, written symbolically as a poly- 
nomial in J, K ; and M turned out to be too complex to be contained in its entirety 
in the system. STEPS 30.9 and 30.95 calculate in turn the coefficients of the 
powers of J from The 0th through the 6th, and these coefficients are listed in 
the array TSTJ (I), 1=0 to 6 as output: 


(A. 16) 
(A. 17) 

(A. 18) 


A-4 



1AM 

WELCOME TO 1AM (72321) 

♦LOAD FROM "PROOF" 

♦ 

DO PART 30 
TSTJC0) i 

8 + K t 9 + K 1 8* ( 48*M + 40') + K»7+(120 *Mt2 + 232 + M + 64) 

+ K t 6* ( i 60 tM t 3 + 560+Mf 2 + 368 + M + 32) 

+ K t 5+ ( 1 20 + M t 4 + 720 + Mt 3 + 872+M»2 + 200 + M - 8) 

+ K*4*( 48+Mt 5 + S20+MT4 + 1088+Mt3 + 48B+MT2 - ) 6*M - 8) 

♦ Kt3+<8*MT6 + 200 + M 1 5 + 752*Mt4 + 600+M»3 + 8*M » 2 - 32*M> 

+ K 1 2+ ( 32*M t 6 + 2 72 +M t 5 + 392+MT4 + 40 + Mt3 - 48*Mt2) 

+ K ♦ < 40 + M t 6 + 1 28 + M 1 5 + 32 + M»4 - 32+MT3) + 16+M*6 + 8 + M»5 - 8+MT4 

TSTJ(l): 

48*K 1 8 + K t 7 * ( 264+M + 200> + KT6+(600*Mt2 + 1036+M + 264) 

+ Kt5*(720+Mt3 + 2208*MT2 + 1296+M + 104) 

+ Kt4*(480*Mt4 + 2472+MT3 + 2560+M*2 + 516+M - 24) 

+ K» 3+ ( 1 68*M» 5 + 1 528 + M* 4 + 2584+Mf3 + 952+M*2 - 4R + M - 16) 

+ K »2+ (24 + Mt 6 + 492 + M t 5 + 1392+MT4 + 816+MT3 - 28*Mt2 - 40 + M) 

+ K * ( 6 4 *M » 6 + 376+Mt 5 + 320+MT4 - 8 + Mt 3 - 32+M»2) 

+ 40+M 1 6 + 44+Mt 5 - 4 + M*4 - 8*Mt3 
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TST J <2 ) 


120*Kt 7 + K » 6* C 588*M + 40 8) + K t 5* ( 1 1 76*M* 2 + 1816+M + 432) 
+ Kt4*( 1224*M*3 + 3248*Mf2 + 1724*M + 128) 

+ Kt3*(696*Mt4 + 2960*M»3 + 2638*Mt2 + 472*M - 24) 

+ K 1 2 * ( 204*M 1 5 + 1432*M»4 + 1908*Mt3 + 588*Mt2 - 40*M - 8) 

+ K*(24*Mt6 + 344*Mf 5 + 638*Mt4 + 280*M»3 - 28*M»2 - 8*M) 

+ 32 *M 1 6 + 7 6*Mt 5 + 36*Mt 4 - l 0*Mt 3 - 2*M»2 

TSTJC3): 

1 60 *K t 6 + K t 5* < 672 + M + 432) + K » 4* C 1 1 28*M* 2 + 156B*M + 352) 
+ K t 3 * ( 960*M 1 3 + 2 1 88 *M t 2 + I 064*M + 72) 

+ K » 2 * C 432 *M t 4 + 1452*Mt3 + 1128*Mr2 + 172 + M - 8) 

+ K*(96*Mf5 + 452 + M t 4 + 474*Mf3 + 108*M»2 - 8*M) 

+ 8*M T 6 + 52 *M * 5 + 58*Mt 4 + 1 1 *M r 3 - 3*M*2 

TSTJC4) : 

1 20*K t 5 + K t 4* ( 408*M + 248) + KT3*(534 *Mt2 + 680*M + 144) 

+ K » 2 * ( 336 *M f 3 + 660*Mt2 + 292*M + 16) 

+ K* ( 1 02*M 1 4 + 260*Mt3 + 170*Mf2 + 16*M) 

+ 1 2 *M t 5 + 32 *M t 4 + 21*M*3 + Mt2 

TST J ( 5 ) : 

48*K t A +- K 1 3 * C 1 20*M + 72) + K 1 2* ( 1 08*Mt 2 + 132*M + 24) 

+ K*(42*M*3 + 72*M»2 + 24*M) + 6*Mt4 + 9*M*3 + 3*M*2 

3 2 2 3 2 

TSTJC6): 8*K + K *(12*M + 8) + K*(6*M + 8*M) + M + M 
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We can conclude that the left hand side of (A. 15) (or (A. 18)) is positive if 
TSTJ(I) > 0 for I = 0 to 6. Recalling that the case m=0 is excluded from this theory, 
it is easily seen by inspection that the coefficients of the various powers of k in 
each TSTJ(I) are positive. Since k > 1 ; it follows that TSTJ(I) > 0 for I = 0 to 6. 

It might be remarked that carrying out this calculation by hand, without 
errors, would indeed be a formidable task - as will be obvious to the reader who 
elected to carry out the omitted proof of the special case (A. 9). 


A -7 



APPENDIX B 


OUTPUT OF WORK1 TO VERIFY CONJECTURE 
FOR m = 3, 4, 5, 6 


_3 

• I AM 

WELCOME TO I AM <7232 1 ) 

♦LOAD FROM "WORKl" 

*M-3i I 1-5>‘Z-5;D0 PART 21 

M: 3 
I 1 : 5 
12: 4 
2: 5 
21: A 

TST1C0): - 77800 
TST 1 ( 1 > : - 132800 
TST 1 (2 ) : 55920 
TST 1 C 3 ) : - 1280 
ODD : 

99328*U t 5 + 3077 1 200*U* A + 379 543 1 424*Ut 3 
+ 2336876 1 3440*IM2 + 7 1 92405 1 47648*U + 88604682485760 

TEST < 0 ) : 

- 77800 + ID5 - 1 9906000*U * 4 - 1 789236000* Ur 3 

- 75 1 35096000*Ut2 - 1 508 764320000* U - 11617066137600 



TEST ( 1 ) 


- 1 32800*1) i 5 - 34068800*1)* 4 - 3238934400*11*3 

- 1 487040 ! 9200*U*2 - 336274364 1 600*U - 30231561830400 

TEST ( 2 ) : 

55920*1) * 5 + 10342560*U*4 

+ 807036480*U*3 + 3263 1 43 1 040*l)t 2 + 670654494720*11 + 5522 714265 600 
TEST C 3 ) : 

- 1 2 80*U * 5 - 49 1 840*U* 4 

- 48885120*11*3 - 2 062344960* Ut 2 - 39 654074880* U - 286123622400 

m - 4 
. I AM 

WELCOME TO I AM (7232 1) 

*LO AD FROM "WORKl” 

*M-4J I 1 -7 ; Z-8JD0 PART 21 

\ 

m: 4 
II: 7 
12 : 6 
Z : 8 
Z 1 : 7 

TST1C0): - 1051359232 
TSTl(l): - 922758144 
TST 1 (2) : 223534080 
TST 1 ( 3 ) : - 15473920 


B-2 



TST 1 C4> : 1 090320 


TST1C5): - 5544 
ODD: 

96680234*1): 7 

+ 1 1 0962352200*U.t 6 + 495733 1 84 ! 0672*1): 5 + 1 1 6557 298 607 1 11 04* U t 4 
+ 15S8276197961984160*U»3 + 1 2 67 7 4 4430258 1 8448000* U t 2 

+ 552404242381 1 804928000*U + 101804347438976286720000 

TEST ( 0 ) : 

- 1 05 1 359232*U» 7 - 669 522022 400*1): 6 

- 1 80675309 10 1 056*Ut 5 - 2 6 869 4 1 1 3 6 60 60032* U : 4 

- 2381 557828659609600*1): 3 - 1 2582555368760 1 1 52000*U»2 

- 36662902 1300957 1 840000*1) - 453 77 077 4 16 69 4087680000 

TEST < 1 ) : 

- 9227581 44*1): 7 - 620726937600*1): 6 

- 1 77442608279552*1): 5 - 280654366462 52544*U: 4 

- 2658741078789734400*1): 3 - 1 5 1 03582 1 06322 1 248000+U: 2 

- 4766493826307727360000*11 - 64479 59 099003068 4160000 

TEST (2 > : 

223534080*1): 7 

♦ 1 3981 5244800*1): 6 + 379 1 2 1 03792640* II: 5 + 57442774098 1 2480*U : 4 
+ 523366306670592000*1): 3 + 2 86 1 384669 7 89 1 840000* 1): 2 
+ 8681 1 5487981 772800000*1) + 11265426294738739200000 


B-3 



TEST C3 ) 


- 1 5473920*1)* 7 - 1 0220672000*11 » 6 

- 27853249 1 7760*U* 5 - 4 1 2389054443520*11* 4 

- 361 1510232 42 24000*11* 3 - 1 8 782 1 06R29 5 1 680000* U* 2 

- 5380068660341 7600000*U - 655112749932748800000 

TEST ( 4 > : 

1090320*11*7 + 560128800*1)* 6 + 1 27579 1 1 9360*U* 5 

+ 1 6490351243520*11* 4 + 1 290 1 02660000000* U* 3 + 60479909582400000 + U *2 
+ 1 562069124000000000*11 + 17068047504076800000 

TEST ( 5 ) : 

- 5544 + U * 7 - 6930000*1)* 6 

- 2 1 1 6743552*U* 5 - 298627426944*1)* 4 - 22862606659200*U*3 

- 98089 1 692704000*U*2 - 22 1 65403863680000*U - 204999081154560000 


m = 5 


I AM 

WELCOME TO I AM (7232 1 ) 
♦LOAD FROM "WORKl" 

*M-5i I 1-95Z-1 liDO PART 21 


M: 5 


II: 9 


12 : 8 


Z: 11 


Z 1 : 10 
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TST 1 ( 0 > : - 17172780156378 
TST 1 C 1 > : - 2 127 1 465952352 
TST 1(2): 3954506831352 
TST 1 ( 3 > t - 284857239744 
TST 1 C 4 ) : 12362967540 

TST 1 (5) : - 290128608 
TST 1 < 6 ) : 8506008 
TST 1 C 7 ) : - 16128 
ODD : 

83277 1 1 488000+U* 9 + 1 5546080348544000+U* 8 

+ 12790995456325632000*IM 7 + 609697046 1 1 09 1 1 3856000*11* 6 
+ 1857536542517692121088000*1)’ 5 ♦ 37 5 45 1 58 4679 455299 1 0 68 1 6000* U* 4 
+ 5038185 88067 3 442 1 93 5 5054 08000* 11* 3 
+ 433070323625221 678285760 1 024000* U* 2 
+ 216481 7394965050586861 51688192000*11 
+ 47967061 6520378954245636463001 6000 

TEST (0 ) : 

- 1 7 1 72 7801 56378*Ui9 - 23484678338339064*1!* 8 

- 1 4 1 96045908490093648*11* 7 - 4975579 169060728536000*11*6 

- 1 1 1 3270555467463550930080 + U* 5 - 1 64 67 9 1 26 1 2 447 42 662082 1 90 08* U * 4 

- 1 607 4 1 0854694929 8820042883 584* U* 3 

- 99564686 1 63043602327 89 67 9 329 28* U* 2 

- 35374636306257631070763277762560*1) 

- 54606409101108981 1207696416768000 
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TEST C 1 > 


- 2 127 1 465952352*11*9 - 3 1 045455262350 1 44*U* 8 

- 20130982577094445056*11*7 - 76 1 387 t 508307544657408*11* 6 

- 1 85 1288387443393074767360*11* 5 - 300 1 1 855705662 1 062644605952*11* 4 

- 324395057 1 7657970 1 3704642 1 S04*Ut 3 

- 2254383345734629 1 1 68784 1 74 1 209 6* U* 2 

- 9 140 179 162 864757607540 658490 5728+U 

- 1 6472 19167649918753516544484966400 

TEST ( 2 ) : 

395450683 1 352*U* 9 + 5636770 1 1 8334432*Ut 8 + 3572809 1 893 1 4805440* U* 7 
+ 132 1041222583083783936*11*6 + 3 1 392006 1 40 1 05533 1 376000*11 1 5 
+ 49708567953334297395666432*11*4 + 52444 10274527423501752707072*11*3 
+ 3 554 4793 8422 681072 1 2378 78 78400*11 1 2 
+ 1 404203 57 9 46812 83349 63289 1 3 551 3 6*U 
+ 246329946269902355895574777036800 

TEST ( 3 > : 

- 284857239744*11 *9 - 398664 1 88733888*11* 8 - 2469 5 5 1 92 9534 60992*1) t 7 

- 88965822613739613696*11*6 - 205 5 1 8563 35936 1 04453 1 20+U t 5 

- 3157832077547419582700544*11*4 - 3227 37 55264098 5564 622 532 60 8* U * 3 

- 21 154547401374726928722051072*11*2 

- 8068491 7 343 89547 095689 7 9009 53 6*U 

- 13640756861286678178141883596800 

TEST < 4 1 : 

12362967 540*11* 9 + 1 634 1 408079920*11* 8 
+ 9637024742785440*11*7 + 33 16981360543804800*11*6 
+ 73274270369561 8862400*11* 5 + 1 0 7 586607508423 1 98749440*11 » 4 
+ 1 048966908774543 1 395947520*11*3 + 6 544625945232 6597 6 7302 75840* U *2 
+ 236980880266 1 7327436482969600*11 + 379288231763912321349058560000 
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TEST (5) : 

- 290 1 28608*1) * 9 - 392029 4 1 2 544*U * 8 

- 22 820961 5682048*11*7 - 7 589 5 1 62363 722240* Ut 6 

- 1 59 79784 1 72 1 07374080*U* 5 - 22 1 5 530833349 5 1 883 4688*U * 4 

- 202563176790072654065664+U* 3 - 1 1 783442334732 50739 5 1 4572 8*U t 2 

- 395775539228071 05281 5564800+U - 5846699438719102726963200000 

TEST C 6 ) : 

8506008*1) * 9 

+ 9142642656*11*8 + 44685284192 64*11*7 + 1 29 1 838248 1 733 1 2*U » 6 
+ 24 131 9265675949440*11* 5 + 29986 1 37704387870208*11* 4 
+ 246471212371721 1088896*11*3 + 1 287025771 73599879962624*U*2 
+ 38632758132809 16627259392*1) + 50686994823099331 195699200 

TEST ( 7 ) : 

- 16128*11*9 - 46928448*U*8 - 32039804 1 60*U* 7 

- 10503523183104*11*6 - 1 99 70 1 43667 7 1 200*U* 5 - 23 632 7888658 6 1 3248*U » 4 

- 1 7687909887524753408*11*3 - 8 1 395564858862428 1 600*U *2 

- 20981 82726 1 430073851 904*U - 231191682528193177190400 

m = 6 

*M-6J I 1<- 13i Z- 15JD0 PART 21 
M: 6 

II: 13 

12: 12 

Z: 15 

Z1 : 14 
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TST l (0) : 
TST 1 Cl): 
TST 1 (2) : 
TST 1 (3) : 
TST 1 (4) : 
TST 1 (5) : 
TST 1 (6) : 
TST 1 (7 ) : 
TST 1 (8) : 
TST 1 (9) : 
TST 1 < 10) 


- 56210426294907000138750 

- 51641842621500873495000 
86831456740844941 12500 

- 578870111235854295000 
22089602975956980750 

- 558494058584919600 
10167863351045400 

- 139563590166000 
1 549199029950 

- 12385933560 

: 122949684 


TST1C11): - 160056 
ODD: 

47160651 19064668823552*Ut 13 + 2906403097 4986852042 7847 68* Ut 1 2 
+ 7691 4586859559054334878449 664*Ut 11 
+ 1 1 8975320098808559566496877969408*Ut 10 
+ 1217472441251 1 5495 1 754 1 47 1 09 07 30 5984* Ut 9 
+ 87763205251204823479 1 660377449 70629 1 20*Ut 8 
+ 4609574625190831 116853019751 48 6550337 1 264* Ut 7 
+ 1 7 92 50367 3658 59 42281819 1 8 680 68 8833 888 727 8592* Ut 6 
+ 51738551778330 4044624686535939770825839529 1 648*Ut 5 
+ 109682799260591431318051 40 7238 79 4567 6009 9 8 60 64384* Ut 4 
+ 1 662306562 5667 00557 499 680282 75 49 6061 448042 483548 1 60*U» 3 
+ 1 7074260918000089340766801 4229 9 550 10829823 569297 40800* Ut 2 
+ 10662446256049939517096498778321 1 4256625735 1 6056 1 664000* U 
+ 30594590881074875512572070686099795417454363772190720000 


B-8 



TEST (0 ) 


- 562 1 042629 4907000 1 38 750*U* 13 - 220746525969335467402372500+Ut 1 2 

- 3996027953 1 4544695887653835000+U* 11 

- 441 4534357 1293 1353 10 10660749 70000+Ut 10 

- 332014522 1 59078 5029 1 88424 1 2 1 8 7300000*11 1 9 

- 1 79492948298305379000218562641 906040000* Ut 8 

- 7 1 749592490147209417478299 1 3 1 2 630 6 7920000* U» 7 

- 214678197291 528 15797 537427 9660 14298 1291 20000*U* 6 

- 48068475772 1 81 5725086661 3654598253985546240000+Ut 5 

- 7953060469680692276719 1903 1 8000329 52 40 7 564800000* Ut 4 

- 944823 1 581 68345694 17 7250379463 5 1092663443845 12 0000* Ut 3 

- 762952 1341287663782655844929 1337208 587 448449 43360000+Ut 2 

- 3751891542831 1111 1 5708 69 532697 7 1 44207967 1 28820 1 2 1 60000*U 

- 84817559855131 14910006183194681202778692427207147520000 

TEST ( 1 ) • 

- 51641 84262 1 500873495000*Ut 13 - 2 1 02 5 1 855 1 4740 558 6499 530000* Ut 1 2 

- 39509294408370623933 6226300000*11* 1 1 

- 45372463967862 1 56557888 1053640000*Ut 10 

- 35527733493297127333016844541 3200000* Ut 9 

- 20030351 1848556393026592559797855840000*U*8 

- 83654480522660340415 111811 00 1 857252 1 60000* Ut 7 

- 2620380460301270681 1 1 03808 1 8 1 3 1 6 1 34460800000+U* 6 

- 615621 60 72 93 12 463509 1232 9 69 6 67141 1 9 69 505280000 + U * 5 

- 1071 37744075697 1375965412321 01 61 7320677893 1200000+U* 4 

- 1342506268700998299948610771 4798523481 5552 798 720000* U* 3 

- 1 1470204007227641899185448593243324499553997946880000*Ut2 

- 5989054188794779761 1 679 4 1 9 7 604847459 623099 039 7 44000000+U 

- 14433277820600385587234154940034794127874994093424640000 


B-9 



TEST (2 ) 


86831456740844941 12500%Ut 13 + 349 60360 7 1 9 97 45460 1 973 5000* Lit 1 2 
+ 6495956328 1965599467355850000*l)t 1 1 
+ 7375 4542153473 881 1 222839 1 5340000* Ut 10 
+ 57089 875 49 7 5 8698 7 69 306 4 5 73 4 60600000% Ut 9 
+ 3181380900724757S974487218328074000000*Ut8 
+ 131307581 58385031400656234939394647520000+UT 7 
+ 4064189807 109491 645087640542 62802 6088 000000% Ut 6 
+ 943334095733240173855032747432925951 188480000%Ut5 
+ 1 621 685373595285406235008737 1 797861 7259801 600000%Ut 4 
+ 20069681 4201 6638963477 1 747 1 508 3 1 05697 527 685 1 200000%U t 3 
+ 1 69324622 1 722 556244 70 6 40429 68 1 2 7304037 476063 64 1 60000%U t 2 
+ 87287851 1 50500630 1 203 1 20 538 1 3 6 1 8467003224843878400000%U 
+ 2076470493414656502362766638956276795005175741808640000 

TEST (3 ) : 

- 5788701 1 1235854295000%Ut 13 - 2302 602402 542 752 1 88450000% Ut 1 2 

- 42 25947780992 1 9 1 868079 1 80000%U T 1 1 

- 4738096959442496902095 129000000%Ut 10 

- 36207740285640778780205563 14000000%Ut9 

- 199147451 5356069 74666976 19 42230240000%Ut 8 

- 81 10574875 1893999 62029 708 7201 2443200000% Ut 7 

- 2 476399872 1 7 52 7 1 62 3 6 59 9 60 7 32 39 62 9 477 7 600 00% Ut 6 

- 5668596272938783323 1 1 04 1 9 6 52 744039 38 56000000%U ♦ 5 

- 960757840258652851 03482937 181278051 484 1 6000000 %Ut 4 

- 1 1 7 1 9120820287500859 1 4998542882670330949304320000%Ut3 

- 974198240421605174490877580855035297349468 1 6000000% Ut 2 

- 494668883061 4721 68194510423635359201 43 5302 49 42080000%U 

- 1 1587141 1764609165297427780196799580236960235520000000 


B-10 



TEST ( 4 ) 


22089602975956980750*11* 13 

+ 86491002543982696588500*11* 12 + 1 56 1 9 5689 4 1 52 1 3 1 30744347000*U r 1 1 
+ 172259143093717771535408250000*11* 10 
+ 12943359645038506261912531 62 60000* Ut 9 
+ 69970518596628 6281707809 19 761 144000*11*8 
+ 2 7 9967 705042407 677 1 9968332 553302032000*11* 7 
+ 839475360731 1 48555 59369 6330 1 0 1 7 40 566 4000*U* 6 
+ 1 8862 62462 178359470453628348 1 062 5 5603200000* U* 5 
+ 3 1367 6 48692 132063429 59 4392 85 63 52 35 76 5785 60000* U* 4 
+ 37 52319931481373801 54232643082489005277-88032000*1)* 3 
+ 305757803845217523491 390604807874858492 1 23 545 6000* U* 2 
+ 15210792 8208571 885762889443855329650933273 1 98592000*1) 

+ 3488954088079574088445883089691182162429083648000000 

TEST ( 5 ) : 

- 5584940585849 19600*U* 13 - 2 1 42 6620 1 9 449 632 48 0800* U * 1 2 

- 37892250741 63506961297600*11* 1 1 - 40898 593 788 59 7339 60658320000*11 ♦ 1 0 

- 30057990080776 5399 57 63043 488000*11* 9 

- 15883661 79 42 530430157 30 43 59 1 15 5200*11*8 

- 620860 71 1883269 1780 1666438901 762 5600*11* 7 

- 1 81 74579 46906934867848048629 3 723909 1200*U* 6 

- 39841805042 5048775 1537 4549 58 1 6 8752 1 280000* U* 5 

- 6459517168573458186339920338039023587328000*11*4 

- 752 81820 71 86657559 70550 689 769 89 6849 6 183 705600* U* 3 

- 597201 6426693473433574908786643065250885 1404800*11*2 

- 28901442104040089869423540788843349532932571 13600*11 

- 6443890613269151455082588471 6384821390501478400000 ■■ 


B-ll 



TEST ( 6 ) 


10 167863351045400*11* 13 + 3 79 8 604 1 39 6482 7 79 600* U» 1 2 

+ 65376576000357477055200*Ut 1 1 + 686239 1 805487 1 62354 1 697600+U* 1 0 

+ 4900967043354531 71 3402 8 790 400* Ut 9 
+ 251 45499320697226640444737 862400*Ut 8 
+ 9 53462066132819 50386 4081 1 3 1 4 1 37600+U* 7 
+ 27049966362316756421312025691 13472000*U»6 
+ 57413 1 509 1 5623 484839 99 52829 3 59 879 577 600+ Ut 5 
+ 9003341345090949 45 87 66 1 4349 86 73289 1 7 50 400* lit 4 
+ 10138443021 47 73292 46 579 5 63 453 409 08909 3632000+U* 3 
+ 77626707230249302371 7 1 0460 1 2 6 1 862 5 1 47207 68000+U* 2 
+ 3621894904707963661 88275330406051949909639 1 6800* U 
+ 777666996997867121 179205957383950265692979200000 

TEST ( 7 ) : 

- 1 39563590 1 66000+U* 13 - 5056078 1 485 1 9 1 1 200*Ut 1 2 

- 842 1 65940795966072000+U* 1 I 

- 854085 1 9389642 1 1 69654400*Ur 1 0 - 58842678843057 1 9 1 2028064000*Ut 9 

- 29082078164150875089076181 7600*11*8 

- 10 6073 02 52 7 50 5592 5 46309 1829 17 1 200* lit 7 

- 2 890 5 7 90 9 7 02 50 48 17 709 47 69 128 4428800* U* 6 

- 5884591266676637266881 562326802022400* U* 5 

- 883793210920424262309414701266292736000*11* 4 

- 95169485783506928321046440862841 1 49849 600* U» 3 

- 69572 14922863672669709835661 61 5899 0827 52000* UT2 

- 30942563564752662384288 1 28 1 45 1 06 586283 5404800* U 

- 6322495954953233870393362639809739436851200000 


B-12 



TEST ( 8 ) 


1 549 1 99029950*11* 1 3 + 53408 50 1 3 7 522900*Ut 1 2 

+ 8480508050939067000*l)t 1 1 + 82027 1 98395483021 25200*Ut 1 0 

+ 5387717233171 13661 357 8 400* Ut 9 * 253605 1 659789230 583822993600* U* 8 
+ 87 972 426757027 4 370784097 500800*U* 7 
+ 22761 576 82 1 461 37051982 10539744000*11* 6 
+ 43912655642784039470723912658739200*11*5 
+ 62372 77 8268829 4460 170480 5 7 579 8 69 18 400* U* 4 
+ 633 83 586361003 43351 541 768403833053 18400* U* 3 
+ 4362839281 1 7 1 8 1 793 73374938 1 38204038758400*Ut 2 
+ 1 8227437 1977464236261 49 0682 4 4720 19 599 3 60000* U 
+ 34900931 454200323183347698382998077440000000 

TEST ( 9 ) : 

- 1 2385933 560*U* 13 - 42094444700880*U* 1 2 - 6488644937387 1 840*U* 1 I 

- 60267 1 71904427442240*U* 10 - 37709 1 69470943679463040*U* 9 

- 1 68056069029293750 1 3866240*Ut 8 - 5492 1 29 4633 6049 3402 6663 1 42 40*U t 7 

- 133306237 63 62801 7! 92 22 72 54 63 040* Ut 6 

- 240340626571584167853323108106240*Ut5 

- 317869 41713 62 62 391 1 4 47 9 086 45085 1 840* Ut 4 

- 2997 1 32 1 675930059997 1 1 58 1 49 1 595182080*Ut 3 

- 190732988639646079303914212051 52792 5760* Ut 2 

- 73405309293761 9389080246 12 32939925 50 4000 *U 

- 128987675480883909365489 172549206016000000 


B-13 



TEST ( 1 0 ) 


122949684*U* 1 3 

+ 356079304 1 52*11* 12 + 47944324818 5 616*11*11 + 39 585 1 65 6 1 5882057 6*U * 1 0 
+ 222779 1 68227694460096*11* 9 + 89 9 489 46380023 64783 57 76* U* 8 
+ 267308930651 9 5852735378 1 76*U* 7 + 5906280793909456838894455296*11 * 6 
+ 9684370571 71 3472 70282 79 1 448 57 6* U* 5 
+ 1 1 620032903 8 16082 5 54 70 129 50 6 1401 6* U* 4 
+ 990 400 1 3 72340066 5478 540800 74 18 67 52* U* 3 
+ 567167099339101496140937881396445184*11*2 
+ 19538681 689 022 03 463 492 3 58 60 49 6 449 53600* U 
+ 305494681042359890085465748301414400000 

TEST ( in: 

- 160056*11* 13 - 787582224*11* 12 - 1 3564483508 1 6*U* 1 1 

- 12612520191 17376*11* 10 - 74 1 9 50 69 1 7829909 7 6* U* 9 

- 298016918468024461056*11*8 - 85039302 5540 7693823 64 1 6* U* 7 

- 1 75643 1 8265935066256 423936*11* 6 - 26352422974576496363439 1 4496*11* 5 

- 2 84 19356761 47 905 1 776593639 8336+U* 4 

- 214292945307251341 1 99 549 75293440*1)* 3 

- 106977 67 1 5 8 69 9 65 3 6 63 62 488 4880 6 7 072*U* 2 

- 3166362330675971 78 74343575 60 4 428800* U 

- 419022432774206248141796356915200000 


B-14 



APPENDIX C 


LAM PROGRAM EVAL AND OUTPUT FOR PARTIAL SUMS 

OF THE POLYNOMIAL PARTS OF P 

nm 

EVALUATED AT £ AND £ FOR m = 3, 4, 5, 6 

1 u 


EVAL Listing 


TYPE EVAL 

!5.1:INZ(2>*2;IN7(3>*5JlNZ<4)*-8JlNZ(5)*llJlNZ<6)*15 
15.2 :$DIST- 1 5 

15.3: FOR L*2 TO 6 . ( M-L 5 Z/ I NZ C L > i K-M 5 Z 1 - Z - 1 J S 
DELETE INZ(L))TYPE M.»K*Z»Z1»A(0)*15BC0)*1;4 
FOR 1*1 TO K, (C*M+I 3 B( I)*- - ( Y-C* CC- 1) )*B( I =1 ) S & 

OK - I + 1 i A ( I )-C*(M + C)*A< I - 1 > ) ; DELETE Cl & 

FOR 1*0 TO K, (A1 ( I >*A( I )*Zr <K- I ) J 4 
A2 ( I >-A<T )*Z 1 t (K-I) ); DELETE Ai« 

FOR 1*0 TO K, B<I >*B C I) *Y» (K- I ) 5 & 

SUM X ( n-Al(K)*B(0HAl(K-n*B(n;TYPE SUMX(.l > J & 
DELETE A1 ( K ) > A 1 ( K - 1 ) } & 

SUMINC1 ) *A2 (K)*B (0)+A2(K- 1 )*B C 1 > J TYPE SUMIN< 1 > i & 

DELETE A2 ( K ) > A2 C K — 1 )>B(0)»B( 1 ) ) 4 

FOR 1*2 TO K, (SUMX( I >*SUMX< I - 1 )+Al (K-I)*BCI )i« 

DELETE SUMXC I - 1 ) , A1 (K- I > J TYPE SUMX(I)i4 

SUM I N ( I )*SUMIN< I - 1 ) +A2 (K-I >*B( I >* 4 

DELETE SUMINCI- 1 ),A2(K-I),BCI )JTYPE SUMIN(I>);4 

DELETE SUMX ( K > > SUMI N (K).M#Z»K»Z1 > 


EVAL Output 


I AM 

WELCOME TO I AM (7232 1 ) 
+LOAD FROM "EVAL” 

*D0 PART 15 

M: 2 


K: 2 


Z: 2 


Z 1 : 1 


C-l 



SUMX< \ *• 

SUMlN C l 7 


SUMX(?-> 


5UM1N (2 


8*V 


+ 96*Y 


: \ 6*Y 


+ 48*Y 


t2*Y 


+ 24*Y + 


l7*Y + 


3<7!*Y 


288 
+ 72 


M*. 3 
Kt 3 
Zi 5 


Z1 5 4 


SUMX 


<i>> - 


+ }0800*Y’ 


SUMlN 


( j ^ : 8640*Y 


SUMX 


(2>* 270*Y 


3 , y + }08000* Y 

. 3600* Y 


SUM! 


N<2>: 288*V 


3 . 5 76»V * ^ 120 * V 


SUMX 


<3>J 145 * Y 


3 cn * - 42000* Y + 

+ 4 l 50*Y . 


SUM 


1N C3): 22 4 * Y 


3 OJfV ? ' - 7 680* Y 

+ 3 392* Y 


C-2 


900000 
+ 460B00 



M! 4 


K : A 
Z: 8 
Z 1 : 7 

SUMX< 1 ) : 

SUM I N ( 1 •) 

SUMXC2) : 

SUM I N ( 2 ) 

SUMX ( 3 ) : 

SUM IN <3 > 

SUMXf A) : 

SUM I N C 4 ) 
M: 5 
K : 5 
Z: 11 
Z 1 : 10 


A 3 

- 2 A 1 92*Y + 1 290240* Y 


A 3 

- 1 6 1 28*Y + 1 1 28960*Y 


A 3 2 

1 88 1 6 * Y - 860160+Y + 25R04800*Y 


4 3 2 

16800*Y - 517 440* Y + 19756800 + Y 


4 3 2 

2 432* Y + 647 1 68* Y - 1 8432000*Y + 412876800+Y 


4 3 2 

5824+Y + 492352* Y - 987B400*Y + 276595200+Y 


4 3 2 

6 528*Y + 40960* Y + 13729792 + Y - 309657600*Y + 5780275200 


4 3 2 

8225*Y + 1 37004*Y + R974252*Y - 146<?41200*Y + 3388291200 


C-3 



4 


5 

SUMX Cl): - 3024000*Y + 199584000+Y 

5 4 

SUMINC1): - 2 4 1 9200* Y + 181440000*Y 

5 4 3 

SUMX C2 ) : 2203200*Y - 176774400*Y + 6586272000*Y 

5 4 3 

SUM IN (2 ) : 1 900800*Y - 129600000*Y + 5443200000+Y 

5 4 3 2 

SUMX < 3 ) : - 1 92 600* Y + 129888000+Y - 609230 1 600* Y + 1 6904 7 64 8000* Y 

5 4 3 2 

SUM I N ( 3 ) : 1 00800+Y + 100800000*Y - 4082400000*Y + 1 27008000000* Y 

SUMXC4): 

539450+Y » 5 

- 1 6522000* Y t 4 + 452R279800*Yt 3 - ] 6 1 53 4 4 1 92 00* Y r 2 + 37 1 9048256000* Y 

SUM I N ( 4 ) : 

600R00*Y t 5 

+ 800000 *Y T 4 + 3 1 7 1 600000*Yt 3 - 9 878 4000000* Y » 2 + 2 540 1 60000000* Y 
SUM X C 5 ) : 

378399*Y»5 + 30182790*Y»4 - 707 1 66 1 08*Y» 3 
+ 121481 147304+YT2 - 3644667290880*Y + 73637155468800 

SUM I N ( 5 ) : 

500800*Y t 5 + 2 9 800000* Y * 4 

- 79200000*Y»3 + 76946400000*Y*2 - 2032 1 28000000* Y + 45722880000000 


C-4 



M: 6 


K: 6 
Z : 15 
Z 1 : 14 

6 5 

SUMXC1): - 547 430400*Y + 43 1 101 44000* Y 

6 5 

SUM I N < 1 ) : - 479001600*Y + 4023 6 1 3 4400*Y 

6 5 4 

SUMXC2): 4 1 4849600*Y - 5 1 1 9 3296000* Y + 2263282560000* Y 

6 5 4 

SUM I N ( 2 ) t 35925 1 200*Y - 4 1 9 1 2640000+Y + 1 97 1 570585600* Y 

SUMXC3) : 

6 5 4 3 

- 1 1 9750400* Y + 3 9688 70 4000* Y - 2766234240000+Y + 90 53 1 302400000* Y 

SUM I N ( 3 ) : 

6 5 4 3 

- 7 5398400*Y + 3 1 977792000*Y - 2 1 1 76 1 285 1 200* Y + 7360530 1 862400* Y 

SUMX ( 4 ) : 

80724600*Yt6 - 1 2 4 3 4 79 6000 * Y t 5 

+ 2 1 87 1 02060000*Yt 4 - 11 3 1 64 1 28000000*Yt 3 + 305543 1 456000000* Y » 2 
SUM I N < 4 ) : 

76728960*Yt6 - 7 57 532 1 600* Y t 5 

+ 1641 1 49959680*Y» 4 - 80965832048640+Y » 3 ♦ 23 1 8567008665 600* Y tg 


C-5 



SUMX (5) : 

26049600*Y t 6 + 7 7 9 49 54000* Y t 5 - 72 75 1 2840000*Y t 4 

+ 90989260200000*Yt 3 - 3888 7 309 44000000* Y t 2 + 9 1 662943680000000* Y 
SUM IN ( 5 > : 

38005632*Y t 6 + 6752309760* Yt 5 - 423 1 1 3209344*Y» 4 

+ 63624906 1 96992*Yf 3 - 2 599 60 5 4339 58400* Y » 2 + 649 l 9876242636800* Y 
SUMXC6): 

37440225 + Y t 6 + 2 07 68 60250* Y t 5 + 4360 1 6722500*Y t 4 - 3 1 694605425000* Y» 3 
+ 3172187731500000*Y»2 - 11 8397968920000000*Y + 2520730951200000000 

SUM I N ( 6 ) : 

45535 1 68* Y » 6 + 2 9 72 482688* Y » 5 + 3460 1 3833984 + Y t 4 - 1 74727 1 4244608* Yt 3 

+ 20678681 5596748S+YT2 - 73936525720780800+Y + 1666276823561011200 


C-6 



